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The time-domain cross-correlation of incoherent and random noise recorded by a series of 
passive sensors contains the impulse response of the medium between these sensors. By using 
noise generated by a can of compressed air sprayed on the surface of a plexiglass plate, we 
J-"- ' are able to reconstruct not only the time of flight but the whole waveforms between the 

sensors. From the reconstruction of the direct Ao and So waves, we derive the dispersion 
curves of the flexural waves, thus estimating the mechanical properties of the material without 
■ a conventional electromechanical source. The dense array of receivers employed here allow a 

precise frequency-wavenumber study of flexural waves, along with a thorough evaluation of the 
' rate of convergence of the correlation with respect to the record length, the frequency, and the 

O distance between the receivers. The reconstruction of the actual amplitude and attenuation of 
the impulse response is also addressed in this paper. 

(accepted for publication in J. Acoust. Soc. Am. 2007) 
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Elastic waves at kHz and MHz frequencies are widely used to evaluate the mechanical properties of structures, 
material and tissues. In plates and shell s, the elastic wa ve equation admits specific propagation modes, denoted 
Lamb waves Rover and Dieules aint ( 2000h; IViktorovl ( 19671 )]. which are related to the traction free condition on both 



-5 

\ sides of the medium. Depending on the purpose of the experiment, different measurement configurations have been 
, i proposed. From the dispersion of Lamb waves, obtained from pitch-catch measures repeate d for different ra nges, one 
has access to the velocities of bulk waves and to the thickness of the plate ( see for instance iGao et all (2003)). Other 



pulse-echo (or impact-echo) techniques have been proposed to asse ss the mechanical properties of the plat e, based on 
the simple and multiple reflection of bulk waves within the plate [Krautkramer and Krautkramerj 1 1990f )] or on the 



resonance of high order Lamb modes Clorennec et al. (12007 )1. Som e studies also concern the dynamic evaluation of 



fatigue and/or crack growth (see for instance Ihn and Chang! ( 2004 ); ling and Finkl ( 1996h ) . In this last application 



several impulse responses of the medium are acquired at different dates and are eventually compared to each others 
to monitor the medium. 



> 

All these techniques require the use of controlled sources and receivers. In the following, they are referred 
£T) . to as "active" experimen ts. Another idea has undergo ne a large development after t he seminal experim e nts o f 
IWeaver and Lobkia ( 200l| ) (see for instance the review of IWapenaar and Fokkemal ( 2006f ). IWeaver and Lobkii ( 20061 ) 
, or iLarosd (|2006l )~ By cross-correlating the incoherent noise recorded by two passive sensors, Weaver and Lobkis 
demonstrated that one could reconstruct the impulse response of the medium as if a source was place d at one sensor. 



This n oise correlation technique (also referred to as "passive imaging", or "seismic interferometry" [Schuster et all 



(120041) ]) requires the use of synchronized sensors, and has the advantage of eliminating any controlled source 



In section II of our paper, we compare the dispersion curves of flexural waves (Lamb waves in the low frequency 
regime) in a plexiglass plate obtained by an active (pitch-catch) experiment to the ones obtained in a passive exper- 
iment. During the passive acquisition, we deployed synchronised r eceivers only, and used a rando m noise source: a 
1 mm thick high pressure air jet produced by a can of compressed air [McBride and Hutchison] (|l976l )] . The experiment 
conducted in a plexiglass plate show at the same time: dispersive (Ao) and non dispersive (So) waves, reverberations, 
and absorption. Since the earth's crust has similar properties for seismic waves, we believe that a plexiglass plate at 
kHz frequencies is a good candidate to built small scale seismic analogous experiments. The question of reconstruct- 
ing not only the phase, but also the amplitude of the wave is addressed at the end of part II. This is of interest in 
seismology, where the role of absorption and attenuation in the correlation has not yet been subject to experimental 
investigation. In section III, we analyze the rate of convergence of the correlations to the real impulse response. The 
role of the record length, the central frequency, and the distance between the receivers is investigated. 

II. EXPERIMENTAL SETUP AND DISPERSION CURVES 

To study actively and passively the dispersion of flexural waves, we built a laboratory experiment using a 1.5mx 1.5m 
large, 0.6 cm thick plexiglass plate. The plate is laid on an open steel frame that supports the edges of the plate but 
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FIG. 1. Experimental setup, (a) In the active experiment (pitch-catch configuration), a broadband piezoelectric source S emits 
a chirp that is sensed by a vertical accelerometer R placed at a distance d. (b) In the passive experiment, a turbulent jet 
produced by a can of compressed air generates a white random noise recorded simultaneously at all sensors Ri. The jet is 
randomly sprayed over the area in gray. 



leaves free the upper and lower sides. The traction-free condition is therefore achieved on both horizontal sides. The 
resulting dispersion relation that connects the pulsation lu and the wave-vector k reads: 



e* \ 9 tan(«ty2 + 7)/ ^ > 

where p 2 — ^- — k 2 and q 2 — ^ — fc 2 ; c s (resp. c p ) is the shear (resp. compressional) velocity, and h is the thickness 

of the plate. The parameter 7 equals for symmetric (5*) modes, and 7r/2 for anti-symmetric (A) modes. In the low 
frequency regime, only two modes are solutions: they are labeled A and S - 

For both acquisitions we used broad-band miniature (3 mm radius) accelerometers (ref.^t 4518 from Bruel&Kjaer). 
They show a flat response in the 20 Hz-70 kHz frequency range. We fixed our accelerometers on our plate using a hot 
chemical glue (phenyl-salicylic acid) that solidifies with cooling (below 43°C). 



A. Active experiment 

In the active experiment, a source S (a piezoelectric polymer) is in a corner of the plate, approximately 30 cm 
from each side and is emitting a 3 s chirp s(t) in a linear range of frequencies / from 1 kHz to 60 kHz. A receiver 
R is initially placed at d —1 cm away from the source toward the center of the plate on a straight graduated line. 
After each acquisition, we moved the receiver a centimeter away from the source down the diagonal; we repeated this 
operation 100 times to cover 100 cm of the plate. The dynamic impulse response hd(t) of the plate is reconstructed 
for each distance d by correlating the record r<i{t) by the source chirp: 



h d (t) = r d (t) x s{t) (2) 

As an example, impulse responses obtained for three different distances d — 10,40,80 cm are displayed in Fig. [5] 
The dispersive Aq mode is dominating the record (a), but the non-dispersive Sq wave (b) and reverberations from 
the edges of the plate (c) are also visible. The transit time in the plate is of the order of a few milliseconds. The 
absorption time of the plate that strongly depends on the frequency is of the order of a few hundreds of milliseconds. 
The signal-to noise in the experiment do not allow to record more than a few tens of millisecond. Records are 
therefore dominated by ballistic waves, but waves reverberated from the plate boundaries are also visible. 

From the set of 100 impulse responses hd(t), a spatio-temporal Fourier (f-k) transform was applied. The resulting 
dispersion curves are displayed in Fig. [31 showing the dispersive Aq mode and the weaker non-dispersive So mode. 
Theoretical dispersion curves are numerically obtained after Eq. [T] (dots), they perfectly fit the experimental data 
for c p = 3130 m/s and c s = 1310 m/s. The dispersion curves are widely used to evaluate the mechanical properties 
of the plate, including thickness, presence of flaws... Nevertheless, the source-receiver configuration is sometimes at 
a disadvantage. As mentioned in the introduction, instead of using the conventional source-receiver configuration, it 
is possible to take advantage of random elastic noise to reconstruct the impulse response hd(t) between two sensors. 
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FIG. 2. Source-receiver impulse response hd(t) for different distances d. (a) Dispersive anti-symmetric Aq mode propagating at 
velocities approximately ranging from 150 to 900 m/s. (b) Non-dispersive So mode propagating at ~2400 m/s. (c) Reflexions 
from surrounding edges. 
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FIG. 3. f-k transform of the source-receiver impulse responses. X-axis: dimensionless wavenumber (h is the slab thickness). 
Dots are theoretical solutions of Eq. [1] 



This is performed in the following section. 



B. Passive experiment 

In the passive experiment, we removed the source and placed 16 receivers separated by 7 cm from each other on 
the graduated line. At the former positi on of the source, w e placed the accelerometer Rq that was kept fixed all 
through the experiment. As proposed by ISabra et al. I (|2007n . we used the noise generated by a turbulent flow. As 
a noise source, we employed a dry air blower. Note that, contrary to Sabra's work, our source was easy-to-handle. 
During the experiment, it was randomly moved to cover a 30 cmx30 cm large area located between the corner of the 
plate and the former active source. Note that the precise knowledge of the noise source is not necessary in the passive 
experiment. We also deployed an array of 16 receivers that allow for precise frequency-wave number analysis, and 
worked at higher frequencies, meaning a much thiner spatial resolution. We focused our attention on the A and Sq 
direct arrivals. Since we wa nt to recons t ruct d irect waves, we chose to spray o nly at one end of the array of receivers 
(end- fire lobes described by lRoux~et" al| (120041) or coherent zones described bv lLarosd (|2006h h Nevertheless, spraying 
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FIG. 4. Passively reconstructed impulse responses (linear color bar, normalized amplitude). 



elsewhere gives the same waveforms but requires much more data, thus much longer acquisitions (see section III). 

The noise in the plate was created by spraying continuously for approximately T=10 s. The 16 receivers recorded 
synchronously this 10 s noise sequence, each record is labeled after its distance d from Rq. Then receivers R\ to i?is 
were translated one centimeter down the graduated line and the 10 s acquisition performed again. This operation was 
repeated seven times to cover the 100 cm of the diagonal with a pitch of 1 cm. We ended with a set of 106 ?"d = o..io5(i) 
records. The time domain cross-correlation between the receiver Rq and the other receivers is processed afterward: 



As mentioned in the introduction, this cross-correlation Crf(r) is very similar to the impulse response hd(t) 
obtained in the pitch-catch experiment. Strictly speaking, the impulse response equals the time-derivative of the 
cross-correlation convolved by the source spectrum (which is almost flat here). Nevertheless, the time-derivative 
operation was not performed here since it does not change the spatio-temporal Fourier transform of the data. The 
series of 100 correlations is displayed in Fig. Q] as a time-distance plot. The dispersive A mode is clearly visible, 
including reverberations at the edges. The symmetric So mode is very weak and almost invisible in this figure. In the 
passive experiment, the size of the noise source is «1 mm large, which is quasi-punctual compared to the Aq and So 
wavelength. In that case, and taking into account the large difference in phase velocities, the energy ratio is clearly 
in favour of the Aq mode, the So mode is much weaker (and additionally less converged). In the active experiment 
the source was 6 mm large, thus exciting the So mode more strongly and the Sq mode was more visible in Fig[3J 

Like in the active experiment, we computed the spatio-temporal Fourier transform of the set of 100 traces Cd(r). 
The resulting dispersion curves are plotted in Fig. The Aq mode is clearly reconstructed, the agreement between 
the active, passive, and theoretical dispersion curves is perfect. Ghosts are also visible in this figure, and are due to 
1) small errors in positioning the sensors 2) the presence of reflections from the edges. As to Fig. [5J the So mode is 
not visible. We therefore zoomed into the early part of the correlations (FigJ5]-left) and muted all the data but the 
weak non-dispersive arrival. We then computed again the f-k transform of the data, as plotted in Fig. [BJright. The 
straight line that emerges from the noise is exactly the So dispersion curve, thus demonstrating that the weak arrival 
in Fig. [BJ- left is indeed the So mode. It is important to emphasize that the weakness of the So mode is not a limitation 
of our correlation technique, but is connected to the noise generation. 



C. Reconstruction of the amplitudes 

Since now, the reconstruction of the phase (the arrival time of the wave) of the Green function by correla- 
tion of noise has been widely studied. Feeble attention was paid to the information carried by the amplitude 
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FIG. 5. Spatio-temporal (f-k) Fourier transform applied to the passive data (linear color bar, normalized amplitude). The 
dispersion curve of the Ao mode is perfectly matching the active experiment. So mode is very weak, probably feebly excited 
by the point-like noise source. The dots show the theoretical dispersion curves as obtained in the active experiment. 




FIG. 6. Left: zoom into the early part of the passive data (saturated color bar). The non-dispersive So wave is pointed by the 
arrows. Right: f-k transform applied to the passive data after muting the Ao mode arrivals (set to zero). The So non-dispersive 
mode is now visible. 



of the reconstructed GF. We can report the first experiments of IWeaver and Lobkis 



(l200ll) who noted that 



both the phase and amplitude of the signals were passively reconstructed, and also lLarose et all (|200bt ) who 
used this amplitude information to study weak localisation without a source. In Sec. III.B( we have seen that the 
phase information in the passive experiment perfectly matches the active data, what about the amplitude of the wave? 

To see if the amplitude decay is recovered using correlations, we plot on Fig. the amplitude decay obtained 
with an active experiment (stars). In a homogeneous and lossless 2D plate, one would expect a decrease of the 
amplitude as oc 1/d - 5 . Because of dispersion and the bandwidth used here (10 kHz), the exponent is slightly 
greater (oc l/d o e ). To be more rigorous, absorption adds another exponential decay. From the active experiment s, 
we found that the absorption time was about 260 ms at 60 kHz and 2000 ms at 1 kHz (similar to lSafaeinili et all fjl996h ) . 

On Fig. we also plot the maximum of the amplitude of the correlations versus the distance d between the sensors 
(circles). Active and passive amplitudes are comparable: the geometrical spreading of the wave is well recovered. 
Nevertheless, small discrepancies are visible and are due the imperfect or irregular coupling between the sensors 
and the plate. Estimating the absorption in the plate from the passive data is therefore more speculative with this 
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FIG. 7. Amplitude decay obtained in an active experiment (stars) and by correlation (circles) without additional processing, 
and with whitening (triangles). Attenuation includes geometrical spreading and absorption. 



experimental setup and will be subject to further investigations. 



Nevertheless, in several applications jRoux et~al\ ( 2004 ); [Sabra et all ( 2005al ): [Shapiro et all (|2005[ )]. an operation 
is performed prior to the correlation in order to balance the contribution of all frequencies. The whitening operation, 
performed here in the 1-60 kHz frequency range, reads: 



f d {t) = IFFT 



FFT (r d {t)) 
\FFT(r d (t))\ 



(4) 



then the correlation fo(t) x r d {t) is evaluated. We compared the active amplitude with the amplitude of the 
correlations using whitened data (triangles in Fig. [7]). This time, though a decay in the correlation is visible, we do 
not recover the dec ay of the actual GF. Note that 1-bit processing (correlate the sign of the records, see for instance 
lLarose et all (|2004l )) give similar disappointing results. Since whitening (or 1-bit) does not preserve the amplitude 
of the records, we stress that correlations will unlikely give the actual attenuation. This is of practical interest for 
applications like seismology. To summarize, in order to recover the amplitude by correlating the incoherent noise, we 
suggest that each record be filtered in a narrow frequency band, and then correlated without additional processing. 
The amplitude decay would therefore be estimated for each frequency, and attenuation derived. On the other hand, 
the whitening (or 1-bit) processing is valuable when the only phase information is targeted (tomography for instance). 



III. RATE OF CONVERGENCE OF CORRELATIONS 



Several authors Campillo and Paul ( 20031 ) : ll)erode et all (|2003l )] have experimentally noticed that the correlations 
C d (r) do not only contain the impulse response (the "signal") between the receivers but also show residual fluctua- 
tions that blur the weaker part of the signal. As soon as residual fluctuations are negligible compared to the reference 
impulse response , we consider that th e correlation s have conv e rged. T he convergence of correlations has been theoret- 
ically studied by I Weaver and Lobkid ( 20051 ) and lSabra et al. | (|2005bl) . We propose here to check the validity of their 
predictions with our experimental data. Of course, theoretical works apply to a perfectly diffuse wavefield. Here, the 
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FIG. 8. Signal-to-Noise Ratio of the correlation versus the duration T of the records. The theoretical fit for each distance d is 
of the form as/T with a the fit parameter. 



noise is generated over a delimited area. Nevertheless, as developed in section III.D, reverberations compensate the 
uneven distribution of source, and we believe that previous theories should apply to our experiment. 



A. Convergence with time T 



To evaluate the degree of convergence of the correlations toward the real impulse response between the sensors, we 
split each noise record into N=1000 sub-records. Each sub-record lasts St =10 ms. Correlations are processed for 
each sub-record and each couple of receivers Rq — Rd- We then compute the amplitude of the residual fluctuations: 



<J d (N, t) = 



' W) - (C d (r)Y 
N — 1 



(5) 



where (.) represents an average over N sub-records. The fluctuations were found to vary weakly with time r, we 
therefore average cr<j(T = N8t, r) over r to get a more robust estimation of the fluctuations <7d(T). To evaluate the 
" signal-to- noise ratio" (SNR) of the correlation, the maximum of each correlation is divided by the residual fluctuations 
ad- Experimental results are plotted in Fig. [8] for three dista nces d and fo r the wh o le 1-60 kHz frequen c y rang e. As 
noticed theoretically and experimentally by several authors Sabra et all ( 2005bf k IWeaver and Lobkisl (2005)], the 
SNR is found to increase like: 



SNR = aVT 



(6) 



where a is the fit coefficient. This means that the longer the records, the better the reconstruction. The SNR 
increases as the square root of the amount of "information grain" contained in the records. This amount corresponds 
to the quantity TAf (where A/ is th e frequency b a ndwid t h) and represents t he number of uncorrelated pieces of 
information transported by the waves jDerode et all 

(fl999h :lL 

arose et all (|2004r )] . In our experiment, a few seconds 

was enough to get the direct arrival of the Aq mode. 

Another point visible in Fig. [5] is the decrease of the coefficient a with the distance d between the pair of receivers. 
This means that the passive reconstruc tion of the i mpuls e response is harder when the receivers are far apart. This 
was noticed in seismology for instance Paul et all ( 2005h : IShapiro and Cam pillo (2004)]. A more precise estimation 
of the convergence with the distance d is developed in the following subsection. 



B. Convergence with distance d 

To study the dependence of the convergence rate on the distance d between the passive sensors, each coefficient a 
was evaluated for each SNR curves (as in Fig. [5]) in the 1-60 kHz frequency range. The resulting a(d) is plotted in 
Fig. [5] and is delimited by two theoretical curves of the form: 
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FIG. 9. The fit coefficient a of the SNR is plotted versus the distance d between the two correlated receivers. 
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FIG. 10. The fit coefficient a of the SNR is plotted versus the central frequency / of the records. 



aocl/d^. (7) 

with 0.5 < (3 < 0.65. At first sight, our result is similar to previous theoretical works that predicted a SNR evolution as 
\fljd in a 2D space. Nevertheless, the theory was developed for a non-dispersive medium, where the wave spreading 
factor is \fl/d. In our dispersive plate, the exponent is slightly greater 0.6, see section II. A). The SNR decay seems 
to be driven by the decay of the actual Green function, note that the role of the absorption can not be evaluated in 
our experiment, and will be subject to further investigations. 



C. Convergence with frequency / 

It is now commonly acknowledged that high frequencies are harder to reconstruct by cross correlation than low 
frequencies. We now propose to quantitatively evaluate the role of the central frequency / of the record in the rate 
of convergence of the correlations. To that end, we performed several whitening operations in different frequency 
band (every 5 kHz, ranging from 1 to 60 kHz). Fluctuations where then evaluated for each frequency band and each 
sub-record. The experimental SNR was fitted by a(f)VT, with a the fit parameter plotted in Fig.[Tnj The prediction 
that SNR oc l/v7 is wei l recovered. 
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FIG. 11. Rate of convergence of the correlation (SNR) with the record time T. The air-jet is randomly sprayed over the dark 
gray areas, doted line is the array of receivers. Upper left plot: sources located on the side of the array of receivers. Lower 
right plot: sources located in others places. 



D. Convergence with the noise sources position 



Several theoretical approaches have been proposed to model the reconstruction of the GF between 
sensors . Some invoke a perfect distribution of sources surrounding the medium to image 



passive 

Derode et al\ (|2003f ); 



IWapenaarl (l2004l. Others refer to a perfect diffuse wavefield, obtained either with scatterers or with thermal noise 
Weaver and Lobkisl ( 200lh ] . Our experiment does not correspond to the first case, since the noise is generated in 
limited area. Does it correspond to the diffuse case? Because ballistic waves strongly dominate the records, this is 
not guaranteed and should be checked thoroughly! 



To answer this point, we chose to spray over an area in the side of the array of receivers (see the dark gray area 
in Fig. Hip . If only ballistic waves were present, the direct waves along the receiver line would not be reconstructed. 
This was not the case here: after an integration over a 35 s long record, we could reconstruct the same waveforms as 
in he active experiment. The only difference is the rate of convergence of the correlations: it was found to be much 
slower (see SNR in Fig. [11] in log-log scale). This proves the elastic wavefield to be (at least slightly) diffuse in the 
plate employed in the experiment. To be more precise, the ballistic and the diffusive regime coexist in our plate. 
Right after a noise source occurs, the propagating wave paquet is in the ballistic regime. It turns into the diffusive 
regime after a dozen of reverberation. The first regime dominates the record; the second is weaker but not completely 
negligible. This last experiment is also a clear evidence that even with imperfectly diffusive wavefield, correlations of 
incoherent noise yield the Green function. This is of importance in applications like seismology where the medium is 
not always perfectly diffusive and the noise sources are not perfectly distributed. Figure [11] also demonstrates that 
when the reconstruction of the direct waves is expected, spraying over an area located in the direction of the array 
requires much less data than any other configuration. 



IV. CONCLUSION AND DISCUSSION 

In this article, we presented two experiments using Lamb waves detected by a series of accelerometers fixed on a 
plate. In the first experiment, a conventional source-receiver configuration was employed to construct the dispersion 
curves of the Aq and So flexural waves. In the second experiment, we only used receivers: we recorded the noise 
generated by a can of compressed air whose turbulent stream generates random excitations. Contrarily to an active 
experiment, the advantage here was that we neither needed to know the position of the noise source, nor to employ 
any electronics for the emission. By correlating the noise records, we reconstructed the impulse response between the 
sensors, and recovered the same dispersion curves as in the active experiment (except for the low frequency So mode 
that was not excited by the thin turbulent jet). 

The amplitude decay of the actual Aq mode with distance was successfully retrieved when correlating the raw 
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records, but was not retrieved using either whitened or 1-bit records. This is of importance for applications where 
passive reconstruction of the attenuation of the media is envisioned. To summarize, whitening (or 1-bit) the data is 
helpful to reconstruct the phase of the GF (for imaging applications or f-k analysis), but should not be employed to 
reconstruct the amplitude. 

Using our dense array of receivers, we also carefully studied the fluctuations of the correlations, which are con- 
nected to the rate of convergence of the correlations to the real impulse response. The Signal-to-Noi se Ratio (or 
correlat i on-to-fluctuation rat i o) ob tained in our experiments is in agreement with previous studies [Sabra et all 
( 2005bl ); IWeaver and Lobkii ( 20051 )]. both qualitatively and quantitatively. We showed that the reconstruction of 



the impulse response is better when: 1) we use long record (more data); 2) we employ close receivers; 3) we work at 
low frequency. The resulting SNR curves were best fitted by: 



The parameter B, as defined by IWeaver and Lobkisl ( 20051 ) was found to be [0.25 ± 0.02] (we take c=35 0m/s as 



a mean value for Aq mode), which is much greater than seismic experiments (Shapiro and Campillc] (|2004l )]. This 
is in part due by the very uneven distribution of sources in our experiment, resulting in an anisotropic wave flux 
that accelerates the convergence of the correlation to the direct A$ wave. This parameter B was found to be ten 
times smaller while spraying on one side of the array of receiver. Concerning the dependence with distance d, the 
decay matches the amplitude decay of the actual Green function. We emphasize that it is not a trivial term when 
dispersive waves are considered. A proper prediction for the SNR should account for the dispersion and the absorption. 

Practically speaking, in the 1-60 kHz frequency range, we obtained very good reconstruction of the impulse re- 
sponses (SNR> 5) for distances up to 100 cm with less than ten seconds of noise. This estimation is fast enough 
to be repeated continuously and provides a route for the passive monitoring of structures and materials. At seis - 
mic frequen c ies, co r relations of ambient noise are alrea dy used to monitor active volca n oes |Brenguier et all ()2007l ) ; 
ISabra etall ( 20061 ): ISens-Schonfelder and Weglerl ( 2006t )]. As suggested bv lSabra et all ( 2007f) in a recent paper, we 



foresee similar applications in the field of on board passive structural health monitoring, in noisy environments such 
as aircrafts or ground vehicles, and also on civil engineering structures as bridges and buildings. 
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